A metric theory of gravity with torsion in extra-dimension 
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We consider a theory of gravity with a hidden extra-dimension and metric-dependent torsion. 
A set of physically motivated constraints are imposed on the geometry so that the torsion stays 
confined to the extra-dimension and the extra-dimension stays hidden at the level of four dimensional 
geodesic motion. At the kinematic level, the theory maps on to General Relativity, but the dynamical 
field equations that follow from the action principle deviate markedly from the standard Einstein 
equations. We study static spherically symmetric vacuum solutions and homogeneous-isotropic 
cosmological solutions that emerge from the field equations. In both cases, we find solutions of 
significant physical interest. Most notably, we find positive mass solutions with naked singularity 
that match the well known Schwarzschild solution at large distances but lack an event horizon. In 
the cosmological context, we find oscillatory scenario in contrast to the inevitable singular big bang 
of the standard cosmology. 



I. INTRODUCTION 

Einstein viewed the space-time as a pseudo- 
Riemannian difFerentiable manifold in order to general- 
ize the special-relativistic flat space-time to include grav- 
ity. This was primarily motivated by the fact that the 
local flatness of the manifold structure naturally imple- 
mented his principle of equivalence. The generalization 
came along with the revolutionary idea that the trajec- 
tory of any freely moving test body is simply a geodesic 
in the curved manifold, and that gravity is not a New- 
tonian instantaneous action force, but an effect of the 
curvature of the space-time manifold. 

The basic constituents of the manifold structure are 
the metric, that defines the distance between any two 
points of the manifold, and the connection, that defines 
the covariant derivative and the curvature of the mani- 
fold. Any theory of gravity should couple the dynamics 
of these quantities to the dynamics of the matter moving 
in the space-time manifold. Among the existing theories, 
the ensuing field equations of general relativity (GR) are 
perhaps the simplest. Since in GR, torsion, the antisym- 
metric combination of connection coefficients, is identi- 
cally zero, and since GR has withstood numerous pre- 
cise experimental tests [I] |2] , introduction of torsion has 
seemed superfluous except in the presence of matter with 
intrinsic spin as in Einstein-Cartan formulations ^3-7 . 

However, for two major reasons, alternate theories of 
gravity that reduce to GR in the weak fleld limit are se- 
riously pursued. The first reason is that GR leads to 
inevitable singularities - black holes (death of a massive 
star) and big bang (birth of the universe). Though it 
is conventionally assumed that quantization would elim- 
inate these singularities, GR is not readily amenable to 
quantization. The second reason is that the standard 
model of cosmology based on GR requires most of the 
universe to be composed of unknown dark energy in or- 
der to account for various cosmological observations [8]. 
A common strategy to construct modified theories of 
gravity is to make the Lagrangian density a nontrivial 



function of the Ricci scalar [51 [TD] and use the action 
principle to derive the modified field equations. Another 
common strategy is to introduce extra-dimensions while 
constraining the physical particles to a (3-f 1) dimensional 
hyper-surface as in the brane- world theories pTHT5] . In 
this paper we explore a different approach by introducing 
metric-dependent torsion in Kaluza-Klein type theories 
[141 115| with one extra dimension. 

In our approach, we consider a five-dimensional (5D) 
manifold foliated by a family of 4D hyper-surfaces, whose 
geometries are virtually indistinguishable from that of 
the 4D space-time of GR. The axis of foliation is spe- 
cial in the sense that there could exist non-vanishing 
torsion components along that dimension. We impose 
constraints on the connection so that any motion in the 
fifth dimension does not affect observations based on the 
geodesic motions along the 4D hyper-surfaces, thus keep- 
ing the fifth dimension essentially hidden. The imposed 
constraints determine uniquely all the non- vanishing tor- 
sion components in terms of the 5D metric fields, making 
this a purely metric theory of gravity. Besides uniquely 
determining the torsion in the 5D geometry, the imposed 
constraints lead to interesting equivalence between the 
5D geometry with torsion and the torsion-free 4D geom- 
etry of GR. In particular, the components of the connec- 
tion and the Ricci tensor along the 4D hypersurfaces turn 
out to exactly match what would arise from GR on a 4D 
space-time. Consequently, any test of this theory based 
on geodesic motions will yield the same results as GR. 

Though by construction, the extra-dimension is hid- 
den at the level of geodesic motion, its effect is clearly 
reflected in the field equations. The field equations are 
obtained by imposing the constraints on the action and 
varying it with respect to the metric. This leads to global 
solutions qualitatively distinct from those obtained from 
GR. Most notably, we find positive mass naked singu- 
larity solutions that match the Schwarzschild solution at 
large distances but lack an event horizon. In the cosmo- 
logical context, we find oscillatory solutions in contrast to 
the inevitable singular big bang in standard cosmology. 



We begin section II with a review of the general frame- 
work of the 5D geometry. Section III deals with the 
specification of the constraints and the determination of 
the torsion and connection in terms of the metric. Sec- 
tion IV is devoted to the derivation of modified Einstein 
equations from the standard action principle using Ricci 
scalar as the Lagrangian density. In section V, we ap- 
ply the modified Einstein equations to the homogeneous 
and isotropic cosmology and identify numerical solutions 
pointing to accelerating and oscillatory solutions to the 
universe. In section VI, we discuss static spherically sym- 
metric vacuum solutions and demonstrate the existence 
of positive mass naked singularity solutions. The final 
section is devoted to a summary and discussion of the 
results. 



II. GENERAL FRAMEWORK OF 5D 
GEOMETRY 

We denote the coordinates of the 5D manifold by the 
latin indices, i,j, k, ... that take values 0,1,2,3 and 5, and 
the coordinates along the 4D hypersurfaces by the greek 
indices, fi, v, A, ... that take values 0,1,2 and 3. Fig. [l]is 
a schematic representation of the 5D geometry. With x^ 
denoting the axis of foliation, the metric of the foliated 
5D geometry has the form: 
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Here A^ is a 4D vector, whose indices are raised and 
lowered with respect to the 4D metric g^" and g^i/. The 
fifth dimension is space-like if e = +1 and it is time-like 
if e = —1. Note that the 5D metric is denoted by the 
bold face g and the 4D metric is light faced g. 

Let us denote the connection in the 5D geometry by 
r. and its antisymmetric part, the torsion by T. . 

T. jfe = '^■jk ~ r. kj- (4) 

Denoting the covariant derivative induced by the con- 
nection by V, the metricity condition is expressed as 
V/cgij = 0. With the metricity condition, the connec- 
tion r. can be expressed as a sum of the Levi-Civita 
connection T and the contorsion K. , 



r' 



jk 



■ ok 



K' 



jfc' 
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FIG. 1. Schematic representation of the 5D geometry. 



where the Levi-Civita connection is expressed purely 
in terms of the metric 

r!,fe = {/'fc } = ^g™ [a^gfe™ + afcg,™ - a„,g,fe] , (6) 



and the contorsion in terms of the torsion [16] . 
1 



K' 



jk 



2 ^^'■'■^ 



t: 



j-fe 



Tfc'.j] 



(7) 



In the absence of torsion, the connection is simply the 
Levi-Civita part. In order to compare the dynamics of 
this geometry to GR, we consider a reference space-time 
in four dimensions with the metric g^^ and torsion-free 
4D Levi-Civita connection r:** . 



r:^^^ = l^g^'^ldf.gay + S^g^a " S^g^z. 



(8) 



This connection is different from the 4D components 
of the 5D Levi-Civit a co nnection that contains additional 
terms f :^ (see eq. |A7| in the appendix) that depend on 
the extra-dimensional metric fields A^ and <I>. Hence the 
4D components of the 5D Levi-Civita connection can be 
written as 



■ ^u 



■ ^u 






(9) 



In the presence of torsion, with the inclusion of contor- 
sion, the 4D components of the 5D connection takes the 
form 



T^A T^A 



r 



/if ' ^^- flw \ ^) 



K' 



We note that the additional terms (f :^ -I- K;** ) do not 
generally vanish. However, in the next section we impose 
constraints on the connection and find that these terms 
do vanish. 



III. CONSTRAINTS ON THE CONNECTION 

With minimal modifications to standard GR in mind, 
we first assume that the 4D components of the connec- 
tion T^ are symmetric, that is (i) T;^ = 0. Next, we 
require that geodesic motion and its observable effects in 
4D are not affected by any motion in the fifth dimension. 



This requirement essentially ensures that the fifth dimen- 
sion stays hidden at the level of 4D geodesies. For this 
purpose, considering the 4D components of the geodesic 
equations in the 5D geometry, namely. 



-f\.x^x'^ + (t\, + r',^) 



x^'x'^ + t%^{i^f 



0, 



we are led to the second constraint (ii) T'^^^ = T^c,^ = 0. 
An alternative formulation of these constraints in terms 
of vielbeins is worked out in [IT] . These constraints are 
clearly not tensorial in nature because the fifth dimension 
is singled out. It turns out however, that they are suffi- 
cient to determine uniquely all the non- vanishing torsion 
components in terms of the metric (see Appendix A for 
details). 
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(11) 



where J^ = ^~^d^/dx\ 

Using the above results for torsion and equations [5j [6] 
and [7] we find the connection coefficients, 
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Here V^ is the covariant derivative operator in the tor- 
sion free 4D geometry with metric g^^. This connection 
has a very special property: its 4D components are ex- 
actly the 4D Levi-Civita connection. That is, F;'' and 



^^ m eq ^ 
In addition 



10 



exactly cancel each other. 

to determining the torsion and the con- 
nection in terms of the metric fields, the constraints also 
imply that the 4D metric on all the hypersurfaces are 
identical. As a consequence, the 4D components of the 
connection also do not depend on x^. 



dg, 



fiiy 



dx^ 



= 



dV: 



fiiy 



dx^ 



= 0. 



(13) 



This should be contrasted with the Kaluza-Klein type 
theories where it is a priori assumed that g^v, A^ and 
$ are independent of x^ , known as the cylindrical con- 
dition. In our framework, though g^^ is required to be 



independent of x 

^5 



A^ and <& can in principle depend on 



Substituting the connection (eq. 12 ) in the Ricci tensor 
defined by 



R^k = dkV- 



■ 3^ 



d,f^. 



X. 



Jl 



-pj -plTl 

• jm ■ kiJ 



(14) 



we find 



^llU ^^U7 



Rij.F, — Rr-,!! — -R55 — 0. 



(15) 

Here R^i^ represents the Ricci tensor constructed from 
the torsion-free 4D Levi-Civita connection. Hence the 
4D components of the Ricci tensor exactly match the 
Ricci tensor in GR with the metric g^j^. It also follows 
that the 5D Ricci scalar is exactly the same as the Ricci 
scalar in the torsion-free 4D space-time, that is R = R. 

An important point to emphasize is that at the level 
of geometry, this framework is virtually indistinguishable 
from the torsion-free 4D space-time of GR. Any observ- 
able geodesic motion or geodesic deviations between par- 
ticles would match what we expect based on GR. How- 
ever, this is true only to the extent the metric g^i/ is 
identical to the solution of the Einstein's equations in 
GR. In section VI, we will see that this is indeed true 
in the weak field limit for spherically symmetric vacuum 
solutions that are relevant for experimental observations 
within the solar system. 



IV. ACTION PRINCIPLE AND MODIFIED 
EINSTEIN EQUATIONS 

We start with the standard Einstein-Hilbert action 
with Ricci scalar as the Lagrangian density. 



S = 



(16) 



In varying the action, we note that the Ricci scalar and 
the connection coefficients described in the previous sec- 
tion are functions of the metric components alone. 



SS 



RS^^ + R,k6g'''V^ 
5R^k g''\/=g d^'x. 



(Fx 



+ /",5i?,fcg'"V=id^a;. (17) 

The first term gives rise to the usual Einstein tensor, 

G.fe = R^k - (l/2)g.fci?. 

In the absence of torsion, the second term becomes a 
boundary integral which vanishes when the variation is 
fixed at the boundary and hence will not contribute to 
the equations of motion. But in the presence of torsion, 
the second term gives a nonzero contribution. 



From eq. 14 we find the variations of the Ricci tensor to be 



5R,, = duStr^ - d,5f\^ + t\jt"], + f™,<5f ^,„ - f ^.^-Jf-, - f ™,<5f ^^,„ 



v,5tr,^ V,<5f^,, 



-L-fcj"-'- • mi- 



Then, the second term in the r.h.s of eq. [17] takes the form 



5R^k^'''V^(Px 



Vfe(g^'=<5f^J - V,{^^St\;) V^cfx + / g^'=T"^/f^"_V^d^ 



(18) 



(19) 



In deriving the above equation, we have used the metricity condition, namely Vjg*'^ = 0. Substituting for the covariant 



derivative in the first term of the r.h.s of eq. 19 



Vkis^'Sfr^) V,(g^^5f\J 



g'fejf:' 
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-gd^x- / g^'St^ 



c'jV=g 



• rri] 



(fx (20) 

^/^cfx. 



The first term in the r.h.s of the above equation is a boundary term, an integral of a total divergence. This will vanish 
when the variation is fixed at the boundary, and hence can be ignored. The second and third terms in the r.h.s can 
be simplified by noting f™„j, = f "J„j. + T"J„j., and f"J„j, = (9fc^/^)/y^, leading to 

Vfc(g''=<5f^",,) - V,{g^Hv\^)] V^d'x ^-J Vl^g^'^Str^^^d':, + J Tr;,„g*'=5f^;,y=irf^x, (21) 

and eq. [19] becomes 



SR^kg"'V^d'^X 



rpm cfikjcf^j I rpm „ik <:f^j , rpm „ifcx-pj 



'gd X. 



(22) 



The second term in the r.h.s of eq 17 is thus a func- 
tion of torsion and variations in the connection given by 
eq. [22] Clearly, this term vanishes if we assume that tor- 
sion is zero, and action principle would yield the standard 
Einstein equations. Alternatively, if we treat the varia- 
tion in the connection to be composed of independent 
variations in metric and torsion, we obtain the Einstein- 
Cartan equations [3], which ultimately leads to zero tor- 
sion when the matter is not coupled to the connection. 
However, since torsion is not an independent degree of 
freedom in our framework and is a function of the met- 
ric components given by eq. |11[ we first substitute its 
components in terms of the metric and then carry out 
the variation with respect to the metric. To this end, we 
note, 

T!i.„g''=5f^,, = T^jg^''(5f«„^ + sr%^) + g'-'sr^.l 



T"I.,g^'=<5r^™ = Tl,[g'^'^sr%^ + g'^'-sr^,] 



(23) 



Taking these terms together, eq. 22 takes the form 



/ 



T^.5 [g^^-^r!-^, - g'-'ST'^^^] V^d'x, (24) 



Note that the variations in the connection, ST^ and 
ST"^ , involve only the 4D components. Since these are 
the 4D Levi-Civita components that only depend on the 
4D metric g^^, the above equation takes the form 



SRik g 



ik 



gd^x= R^Jg^"" ^^d'x, (25) 



where (see Appendix B for details) 



M-v - V(^B,) - (V • B)g,,, + J(,,B,) - (J • B)g^„ 
B„. = T^,, = J„ - d.A., - A, J 



(26) 



Taking together the variations in both terms in eq. |17[ 
we obtain the modified Einstein tensor. 



tfpiy — R 

G55 



^i^ 



n ISa"' 



A„A^e$2)^_^jj,,_, = E.^, 



ie$2^ = I]55, 



where E is the stress energy tensor that arises from the 
variations of the assumed matter fields in the Lagrangian. 
Our focus being on gravity, we will not further discuss 
the origin of E. Since the physical interpretation of the 
stress- energy is more transparent with one covariant and 



one contravariant indices, we express the above equations 
in an alternate form, by noting G/ = g'^^Gik, 



HOMOGENEOUS-ISOTROPIC 
COSMOLOGY 
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^5 ' 
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(27) 
(28) 
(29) 



These are the modified Einstein equations in our frame- 
work. Since by construction the fifth dimension is hidden 
with respect to the observable 4D motion, the 5D compo- 



nents of the stress tensor S.f 



and Yi^ are unobservable. 



It is not possible to solve equations [28| and [29] unless these 
components are theoretically known from the 5D matter 
Lagrangian. In the present formulation, for simplicity, 
we shall ignore these equations as though they simply 
serve to evaluate the components I]„^ and £5^, and treat 
only eq. [27] with the observable 4D stress tensor to be 
relevant to physical solutions. In the absence of specified 
matter fields in the Lagrangian, an alternate way to in- 
terpret the modified Einstein equations is to regard — H " 
as extra- dimensionally induced matter. 

When H,," = 0, eq.[27 reduces to the standard Einstein 

In this case 



-M - 0, eq ^^ 
equations for the 4D metric components g^i. 



the 4D Bianchi identity necessarily implies the conserva- 
tion of matter VjyS/' — 0. But in general when H^ is 
non-vanishing and dependent on the extra-dimensional 
metric fields A,, 



and <I>, eq. 



27 by itself may not be suf- 
and <I>, even after 



ficient to solve for g^^ along with A^ 
fixing the gauge. However, an important physical simpli- 
fication can be achieved by generalizing the cylindrical 
condition to assume that A^ and <& do not depend on x^. 
With this assumption, B^ — J^ and H "^ depends only 



on $ and not on A^, and eq. 27 is sufficient to solve for 
both g^j^ and $ 
uated from eq. 



28 



The vector A^ can in principle be eval- 
by setting S^ to zero, but this would 
be inconsequential as A^ is decoupled from the physi- 
cally relevant equation that solves for the 4D metric g^,y. 
Hence, in the rest of the paper we will make the assump- 
tion of cylindrical condition in order to explore solutions 
of physical interest to the modified Einstein equations. 

Finally, when H " is non-vanishing, we note that S '^ 
does not necessarily have to satisfy the 4D matter con- 
servation. However, with minimum modifications to GR 
and the empirical conservation laws in mind, it is reason- 
able to assert the conservation of E,''. Since the standard 
Einstein tensor satisfies the 4D Bianchi identity indepen- 
dently of H ", the 4D matter conservation implies. 



The 4D metric of a homogeneous and isotropic universe 
has the form 



ds^ 



-dt' 



a^(t) 



dr^ 



1 — fcr^ 



r'^dn'^ 



(31) 



The values of A; = 0, +1,— 1 correspond respectively to 
flat, closed, and hyperbolic spatial geometries. The stan- 
dard Einstien tensor for this metric is given by [19] 



G\ = 3(d/a)2 + 3fc/a2, 

Gl = 2{a/a) + {a/af + k/a^, 

GO r~\(k c^r 
fi — ^A. — ^^-ri 



(32) 



where over-dot denotes a derivative with respect to 
time. Since the geometry is spatially homogeneous and 
isotropic, the metric fields including A^ and $ in the 5D 
geometry only depends on time. Hence the only non- 
vanishing component of J^ is Jt. The induced matter 
terms given in eq. [26[ are 



H, 



3Jt(a/a), 



H/-2Jt(d/a) + Jt + j2, 

Ho TT TT »" 
a — tl^ — n^ • 



(33) 



Before writing out the modified Einstein equations, we 
note that the conservation equation (eq. 30 ) now reduces 
to 



V.H; = 3Jt [{a/a) - {a/a)h] = 







Jt = d/d 



(34) 



These are the only two possibilities. From the defini- 
tion of J^, this implies that either $ is a constant, which 
would give rise to the usual Friedman-Robertson- Walker 
cosmology, or $ = d(i). Focusing on the latter case, H "^ 
simplifies to 



H/ == 3d/o 



H„ 



2(d/a) + (a/d). (35) 



Taking the stress tensor to be that of a perfect fluid, 
the modified Einstein equations (eq. [27|) take the form 



3(d/a)2 + 2,k/a^ + 3d/a = 87rp (36) 

4(d/a) + {a/af + (a/d) + k/a^ = -8ttP (37) 

where p and P are the density and pressure of the 4D 
matter. Combining the above equations, we find 



v.s; = : 



v.h;- = 0. 



(30) 



In the reminder of the paper, we study the solutions to 
the modified Einstein equations (eq. 27 1 in two extremely 
symmetric situations, namely, the homogeneous-isotropic 
geometry and the static spherically symmetric geometry. 



p + 3{p + P)d/a = 0, 



(38) 



which is just a restatement of the 4D matter conservation 
equation. For matter dominated universe, P = 0, and 
consequently eq. [38[ yields 



pa = constant 



Pole? 



(39) 



k=0,q^ = -0.2 



k=0, q^ = +0.4 



k=0, q„ = +0.5 




FIG. 2. Solutions for a{t) for different values of deceleration parameter for k = Q and k — +1. 



In effect, eqns. 36 and 



which needs two initia 
fication of p^- 



39 are sufficient to solve for a{t), 



and a minimum a„ 



The acceleration reaches a 



conditions along with the speci- 



a d + a{a) + ka ^ 87rpo/3, 



(40) 



Without loss of generality we choose the current epoch to 
be at t = 0, set the current size of the universe a(0) = 1, 
and the unit of time such that a(0) = 1. In this unit of 
time, the current value of the Hubble's constant will be 
1. Rather than specifying the value of po, we choose to 
specify the current value of a. The effective equation for 
a{t) then takes the form 



a a + a{d) + ka = 1 + k — Qo, 



(41) 



where Qo = — a(0)a(0)/a^(0), the current value of the 
deceleration parameter, is the only free parameter to be 
specified. Figure [2] shows the behavior of a{t) for vari- 
ous values of Qo for spatially flat and closed topologies. 
For the spatially flat topology, shown in the top panels 
of the flgure, we find that the universe does not origi- 
nate from a singular big bang for all Qo < -1-0.5. For the 
spatially closed topology shown in the bottom panels, we 
find oscillatory solutions for all go < +1. 

Oscillatory solutions in the spatially closed topology 
exhibit a scale factor that oscillates between a maximum 



positive value in a narrow interval around «„„, and then 
becomes negative for the rest of the cycle until it gets 
back near amin- By taking the value of qo arbitrarily 
close to 1, we can make amm arbitrarily close to zero. 
This can be seen from the bottom-left and bottom-middle 
panels of fig. [2] Thus one could construct a universe that 
collapses and bounces back to expand when it reaches 
an arbitrarily small size, or equivalently arbitrarily high 
energy densities. It remains to be seen if such solutions 
would fit the empirical red shift data. 



VI. STATIC SPHERICALLY SYMMETRIC 
VACUUM SOLUTIONS 

The most general static spherically symmetric 4D met- 
ric has the forrrQ 



ds^ = -A{r)dt'^ + B(r)dr'^ + r'^dVl^, (42) 



^ The scalar functions A{r) and B{r) defined in this section should 
not be confused with the vectors A^ and B^ defined in eqns. [s] 
and 1261 respectively. 



and the standard Einstein tensor for this metric is [12] 
rB' + B'^ - B 



GH 



Gr 
r 






^252 
AB -rA' - A 

r^AB ' 
2A^B' - 2ABA' 



2rABA" + rBA''^ + rAA'B' 



ArA^B^ 



(43) 



where a prime in the above equations denotes a derivative 
with respect to r. 

The additional term H "^ in the modified Einstein equa- 



tions (eq. 27) depends only on $ when the cylindrical 



condition is imposed on all metric components and is 
given by 



h; = v^j'' - (V • 3)sr, + j^r - (j • 3)6^,. 



(44) 



Since J^ = <i>^^9p<f>, the quantity V^J'' is intrinsically 
symmetric in fi and ly. The static spherical symmetry of 
the geometry implies that Jr is the only non-vanishing 
component, which we denote by J(r). With this, 

rJB' - 2B (r.r + 2J + r.P) 



H/ 



H_ 



H / = H 



2rB2 
-J{rA' + 4A) 
2rAB ' 
rJAB' - 2rABJ' - rJBA' - 2ABJ (1 + rJ) 



2rAB^ 



(45) 



In order to obtain vacuum solutions, we set Y^f,=0 in 
eq. [27] and find the following three equations. 

J{l + B) 



J'--- 



A' = - 



2A{1 - B + 2rJ) 



B' 



r{2 + Jr) 
2B(rV2 + (l + rJ)(l 



B)) 



(46) 



r{2 + Jr) 

A close examination of the above equations reveals two 
basic properties of the function J{r). First, if J{r) is a 
constant, it has to be identically zero. Secondly, if J{r) 
vanishes at some point, it has to vanish identically every- 
where. The simplest solution to the coupled equations 



(eq. 46 1 is when J{r) vanishes everywhere. 



J(r) = 0,A(r)=(l-^),B(r)=(l-^ 

. ' (^^) 
which of course is the well-known Schwarzschild solution 

as expected. 



A. General solution to J{r) 

Let F{r) = \/rJ{r) when J{r) is non-vanishing. Sub- 
stituting for J{r) in terms of F[r), the coupled equations 




FIG. 3. /3 is plotted as a function of c. 



(eq. 461 lead to the following equation for F{r), 

r [(2F2 -t- F)F" + F'2] = F'{F + 2), (48) 

whose solution in turn determines the 4D metric func- 



tions A{r) and B{r). The obvious solution of eq. 48 
F{r) =constant. This leads to 



IS 



Jir) = c/r , A{r) = {r)-^Tr , B{r) = 0, (49) 

which is clearly unacceptable because B{r) is identically 
zero. 

Assuming that F(r) is not a constant, we can ob- 
tain solutions to the second order differential equation 
(eq. 48 1. In principle, the solution would have two in- 
tegration constants that would be determined by the 
boundary conditions, one of which immediately follows 
from the form of the equation. It can be easily seen that 
if F{r) is a solution, then F{Xr) is also a solution for 
any scaling constant A. We find a general solution in the 
implicit form 



1+/3| 



X'r' = Y \Fh + [3 + ir^ |F/7 + /? - 1| 



:1-I3 



(50) 



where /3 and 7 are defined in terms of an independent 
arbitrary constant c. 
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= Vl + c + c2, 13 = 



1 + c 

Vl + c + c^ 



(51) 



With A and c as two arbitrary constants, eq. [50] repre- 
sents the general solution to the second order differential 
equation (eq. 48 1 . From eq. 50 we find the derivatives of 
F{r) to be 



F' = 



IF^ + 2{l + c)F- 
r F 



c F' 
^" = --4, (52) 



and substituting them in eq. 46 we obtain the metric 



functions A{r) and B{r) in terms of F(r). 

A'{r) = 2cA{r)/rF{r), (53) 

B{r) = l + 2{l + c)/F{r)+c/F^{r). (54) 

In order to obtain asymptotically flat solutions, we 
shall impose the boundary conditions A{r — >■ 00) = 1 




FIG. 4. Numerical solution to F{r) with A = — 1. 



and B{r -^ oo) = 1. To understand the behavior of the 
functions A{r) and B{r) which define the observable 4D 
geometry, we start with the properties of F{r). 

Unfortunately eq. [50] does not yield an explicit func- 
tional form for F(r) except for simple cases when c is 
either or -1. Nevertheless the relevant properties of 
F(r) can be inferred from analyzing this implicit func- 
tion. First note that the quantity (1 + c + c^) is posi- 
tive definite and f3 is finite and bounded for all values 
of c. Figure. [3] plots the behavior of /3 to show that it 
asymptotically reaches +1 and -1 at c = +oo and — cxi 
respectively. The following observations summarize the 
qualitative properties of F{r). 

1) At r = 0, F can take one of two possible values. 
If c < (/3 < 1), then F{0) can be either -f{-(3 - 1) or 
7(-,3 + 1), while if c > (/3 > 1), then F(0) can only be 

7(-/3-l)- _ 

2) In the limit r ^ oo, F necessarily has to diverge in 
ord er t o satisfy the boundary condition B{r — ^ oo) — 1. 
Eq. 50 then implies for large r, X^r^ — F'^, implying F 
could be either positive or negative, such that 



F{r ^ oo) = Xr ^ J{r ^ cx)) = A ^ /r^ 



(55) 



The behavior of F(r) at the extremities is summarized 
in the following table. 





c < (/3 < 1) 


c> (/3 > 1) 


r = 


\7(-/3-l) <0 


F = 7(-/3-l) <0 


7' ^> oo 


F = Ar 


F ^ Xr 



3) From eq. 52 it can be shown that F' — when F 
is either 7(— /3 + 1) or "f{~/3 — 1), which can happen only 
at r = 0. Hence F' is either positive definite or negative 
definite, and so F(r) is a monotonic function. 

4) Though F is monotonic and finite for any finite r, 
it can reach zero at Tq given by 



\2 2 

A r„ = 



/3 + 1 



/3-1 



/9 



(56) 



and from eq. 52 F' diverges at r^. Considering the phys- 
ical relevance of these solutions, we shall only focus on 
solutions that are non-vanishing everywhere. Such so- 
lutions do indeed exist for a range of parameter values. 
Rewriting eq. [50] at r = Tq, 



7(/3 + l) 



i+P 



7(/3-l) 



1-/3 



= 1, 



(57) 



we note that F{ro) = is not the only solution. Nu- 
merical plots in fig. [4] demonstrates the existence of non- 
vanishing F{r) solutions. 

5) From the table above, (i) for c > 0, since F{Q) is 
negative, F{r) has to be negative definite which requires 
A to be negative, (ii) For c < 0, A can be either positive or 
negative, making Fir) either positive definite or negative 
definite respectively. 

The functional form of F{r) described by the above five 
properties along with equations [53] and [54] will yield the 
functional form of the metric functions A{r) and B{r). 



B. Metric functions A{r) and B{r) 

With the boundary condition yl(oo) = B{oo) = 1, 
eqns. [53] and [54| yield 



A{r) — exp I — / 



2c 

rF{r) 



dr 



B{r) = 1 + 2(1 + c)/F{r) + c/F'^{r). 



(58) 



(59) 



The following observations summarize the qualitative be- 
havior of A{r) and B{r). 



c<0, X>Q; M<0 i 
I c>0, \<0; M<0 




B(r) 




FIG. 5. Schematic behavior of A{r) and B{r) 



1) From the asymptotic behavior oi F{r) — > Xr for 
large r, we note that 

9„\-l 

^(r) = l--^^— +0(l/r2), (60) 

r 



Hence, when \c\ 3> 1 and r — > oo, the above solutions 
approximate the Schwarzschild solution with mass M = 
[cA~^]. When both c and A are either positive or negative, 
the gravity is attractive, while when one is positive and 
the other is negative, the gravity is repulsive. 

2) Since F{r) is either positive definite or negative def- 
inite, both A{r) and B{r) are finite and positive for all 
r > 0. At r = 0, since F{0) is either 7(— /? — 1) or 
7(-/3 + l), eq. [59] implies 5(0) =0. 

3) As r -^ 0, the integral in eq. [58] diverges as 
[2c/i^(0)]ln(r). When [c/F{0)] is positive, then A{0) = 
0, and when [c/F(0)] is negative A{0) = oo. The sign 
of [c/i^(0)] is the same as the sign of M ~ [cX^^]. For 
M > 0, A{r) monotonically increases from A{Q) = to 



yl(oo) = 1; for M < 0, A{r) monotonically decreases 
from A{0) = +oo to yl(oo) = 1. 

4) Irrespective of the sign of M, B{0) = and B{oo) — 
1. However, B{r) is not necessarily monotonic. From 
eq. [59] we see that B' = when either i^' = or when 
F{r) — -~c/(l + c). From the previous subsection, F' 7^ 
for all r > 0, but F{r) could attain the value — c/(l + 
c) for certain values of c and A. Since F(r) is a non- 
vanishing monotonic function taking all values from F(0) 
to ±cxi, it is straightforward to check if it would attain 
the value — c/(l + c). When A < 0, F{r) is negative 
definite, and — c/(l + c) needs to be a negative number 
lesser than F (0) = j{— /3 — 1) , which happens only when 
c < —1. When A > 0, F{r) is positive definite, and 
— c/(l + c) needs to be a positive number greater than 
F(0) = 7(— /3 + 1), which happens only when — 1 < c < 0. 

The qualitative behavior of the functions A(r) and 
B(r) for the various allowed ranges of c and A are shown 
in figure^ When M = [cX^^] is positive, A{r) is a mono- 
tonically increasing function leading to attractive grav- 
ity, and corresponds to Schwarzschild solution at large r 
when IcI ^ 1. 



C. Naked singularity at r = 



turns out that 



An important point to note from figure [5] is that these 
solutions do not have an event horizon because both A(r) 
and B{r) arc finite and positive for all r > 0. 

Clearly, these solutions are smooth for all r > 0. How- 
ever, the point r = is a physical singularity. Explicit 
calculations show that the Ricci scalar i?'^ vanishes ev- 
erywhere, but the quantity R^vR^'^ is nonvanishing. It 



RavR'^^ 



2{3F^{r)+2cF{r) 



F\r)B'^{r)r^ 



-c^) 



(62) 



At r = 0, the numerator does not generally vanish, 
but the denominator vanishes, making R^^R^^ diverge. 
Hence these solutions correspond to a naked singularity 
at r = with no event horizon to censor it. 
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VII. SUMMARY AND DISCUSSION 

Metric and torsion are two independent constituents 
of metric compatible Riemannian geometry. Because of 
the immense successes of torsion-free GR, torsion has not 
played a significant role in theories of gravity. However, 
when gravity is to be included with other interactions of 
elementary particles with intrinsic spin, a more general 
theory including torsion becomes imperative [5]- 

In the present work, torsion is incorporated in a novel 
way in higher dimensional Kaluza-Klein type theories. 
Here torsion is not an independent degree of freedom cou- 
pled to spin, rather it is determined in terms of metric 
through a set of physically motivated constraints, which 
serve (i) to confine torsion to the extra dimension, leav- 
ing the 4D space-time torsion free, and (ii) to ensure that 
geodesic motions in 4D remain totally unaffected by the 
presence of the extra-dimension. These constraints have 
previously been imposed in terms of veilbeins |171 I18j . 
but here it is realized that they impose essentially the re- 
quirement that the fifth dimension is hidden at the level 
of geodesic motion. It turns out that the non-vanishing 
torsion components are functions of the 5D metric com- 
ponents with the 4D metric g^j^ obeying the so called 
cylindrical condition, namely, it is independent of x^. 
In the resulting geometry, all the 4D hypersurfaces are 
equivalent, and the 4D components of the connection and 
the Ricci tensor exactly match those of the standard 4D 
GR. Hence, at the level of geodesies, this geometry is vir- 
tually indistinguishable from that of the standard GR. 

To proceed further, we derive modified Einstein equa- 
tions from the action principle with Ricci scalar as the La- 
grangian density. In this respect, an alternate approach 
presents itself. In the Palatini formulation of GR, the 
action is varied independently with respect to the con- 
nection, and in the absence of torsion, the metricity fol- 
lows from the equations of motion. Recently [201 EI], 
it has been shown that even without assuming the ab- 
sence of torsion, variations of the action independently 
with respect to the metric and the connection lead to 
GR equations along with metricity, provided the mat- 
ter Lagrangian is not coupled to the connection. In our 
case, with the constraints imposed on the connection, it 
is more convenient and natural to impose metricity prior 
to action variation. Since the entire connection is de- 
termined to be a function of the metric, we only need 
to vary the action with respect to the metric to obtain 
the modified equations, making the theory a purely met- 
ric theory of gravity. However, adopting a Palatini-style 
approach, one could relax the assumption of metricity 
and vary the action independently with respect to met- 
ric and connection along with the imposed constraints, 
which might lead to a different set of modified Einstein 
equations. 

We apply the ensuing modified Einstein equations to 
study the cosmology of a homogenous-isotropic universe. 
In matter dominated phase of the universe (zero pres- 
sure), we obtain a second order differential equation for 



the scale factor a{t) in contrast to the first order differ- 
ential equation in the usual FRW cosmology. In FRW 
cosmology, the second derivative of a{t) cannot be inde- 
pendently prescribed as an initial condition, and deceler- 
ating expansion is a necessary outcome in the absence of 
a cosmological constant. However in our case, we have 
the choice of an initial condition for the second derivative 
of a{t), which can be tuned to fit the observed acceler- 
ation of the universe. Figure [2] presents the behavior of 
a{t) for various choices of the current acceleration. 

In an earlier version of the present work [22] , the field 
equations were derived differently; by varying the action 
with respect to the metric prior to expressing the Ricci 
tensor in terms of the metric. In the cosmology equa- 
tions generated from those field equations, acceleration 
was not an independently prescribable initial condition. 
Chen and Jing ^3\ showed that those equations yield 
accelerating universe solutions without resorting to dark 
energy. They show that the model not only fits the su- 
pernovae data, but also solves the cosmic age problem of 
old high redshift objects ^^. Whether the cosmological 
solutions described in the current work would fit empiri- 
cal results just as well needs to be investigated. 

In the case of spherically symmetric vacuum solutions 
to the modified Einstein equations, we find some remark- 
ably interesting results. As is well known, in the unique 
Schwarzschild solution of GR, when the mass is positive, 
an event horizon censors the central singularity. In con- 
trast, we find positive mass naked singularity solutions 
that lack an event horizon. Recently, similar positive 
mass solutions without a horizon have been found [25] in 
a simpler setting of torsion-free GR with multiple extra- 
dimensions. It would be interesting to see how the partic- 
ular vacuum solutions in the torsion induced geometry in 
the present case match with those solutions in the torsion 
free geometry. 

The existence of positive mass solutions with naked 
singularity have immediate consequences on gravitational 
collapse, opening up the possibility of an arbitrarily large 
star collapsing to an arbitrarily small non-singular state. 
Since trapped surfaces would not necessarily form in such 
collapses, finite matter pressure could be sufficient to 
withstand a total collapse to singularity. This suggests 
a detailed analysis of such solutions by treating — H "^ in 
eq.[27]as extra-dimensionally induced matter in standard 
GR. It then raises the possibility of a sufficiently strong 
gravitational collapse that stops short of collapsing to a 
singularity with a finite induced stress-energy tensor that 
potentially violates the weak energy condition in the re- 
gion near the center. It remains to be seen if an arbitrar- 
ily small static model star with finite stress-energy tensor 
can be constructed with the external geometry matching 
the type of solutions discussed in this paper. 

In conclusion, inclusion of torsion in the context of 
extra-dimensions presents a novel way of obtaining modi- 
fied Einstein equations that have significant physical con- 
sequences. For clarity and simplicity, we have confined 
the treatment to five dimensions. However, the frame- 
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work can be generalized to arbitrary dimensions D, pro- 
ducing torsion-free D — 1 dimensional metric theory. We 
could also consider many extra-dimensions and general- 
ize the constraints so as to hide all the extra-dimensions 
using torsion, leading to a more general theory. 
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Appendix A: Computing the geometric quantities in the 5D manifold 



In this appendix, expression for torsion, connection coefhcients and the Ricci tensor of the 5D geometry are computed 
in terms of the metric. It will be shown that the physical constraints imposed on the connection will be sufficient to 
determine uniquely all the components of torsion and hence the other geometric quantities in terms of the metric. 



1. 5D Levi-Civita Connections 



The 5D Levi-Civita connection is given by 



^Uk = oS''™ (^fcgim + ^igfem " ^mgifc) 



(Al) 
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Expressing the 5D metric in terms of the 4D metric and the extra-dimensional metric fields given by eq. [3] we find 
the 5D Levi-Civita connection to be 

- ^ (A^Aa + 6$-") ^5 (eA^A,$2) _ 1 (a^Aa + e*-^) Q^g^^ (A2) 

r^r, = -e'^^J" + f^^dr.A" + e^^AMjg _ e'^^A^d^g^'' (A3) 

t%^ = e$2A^ (Ja - ^sAa) - e$2A^AAJ5 + J5 (A4) 

f '^5. = \^^^n ~ e^^A,,r + ig^^95 (e$'A,A,) + ^g^'^agg.. (A5) 

V%, = ^e^'A^FA. + e$'A,A^JA + J. - \A^d^ (eAaA^^^) _ lA^^sg.A (A6) 

^\. - ^e$2 (A,F^^ + A^F,^) - 6$2A^A.J^ + ^A^^s (eA^A.$2) + ^A^^sg^, (A7) 

where F^,^ = d^^Ay — S^A^, J^ = e$^^9ie$ and F;^^^ is the 4D Levi-Civita connection obtained from the metric g^^. 
and V^ is the derivative operator with the 4D Levi-Civita connection. The raising and lowering of indices on F^i,, 
A^ and J^ are performed with respect to the 4D metric. 

2. Torsion components 

We start with the conditions that T:''^^^ = and T^^^ = T^^^ — 0. Together these conditions imply T^^^ — 0. 
The remaining non- vanishing components of torsion that need to be determinned are T^^ -, a total of 10 independent 

components. The conditions T%i = 0, consisting of 20 equations are sufficient to determine uniquely all the non- 
vanishing torsion components. 

f^.5=r'^i5+K'^.5=0 (A8) 

From the above equation, the non vanishing components of contorsion and torsion can be determined in terms of 
components of the components of 5D Levi-Civita connection F. . First taking i = 5, 

K^55 = g^^ §55X^,5 = g'^'^gssT^^.s = -r'^ss (A9) 

From the 5D metric (eq. [2]), we have g^'"' = g^'^ . Multiplying both sides by g^o- and using the orthogonality relations 
of the metric, 

gM^g^'^gSsT^s = -g^.f ^55 

g55T.£,5 = -g^crL^jg 

=> Tl^^s = -g^.f'^sge*-^ (AlO) 



Substituting for f ^55 from eq. 



A3 



we obtain 



T^5 = gM-[J^ - d^A^ - A^Js + A,95g^1 (All) 

We have thus far used 4 equations and solved for 4 of the 10 independent torsion components. Next, take i = u m. 



eq. 



A8 The contorsion components Kf'j^g can be expressed in terms of the torsion components as follows: 



= ^g^' (gssT^s.) + ^g^'^ (gssT!^.. + g.sT^s) 

= i^A^^^T^s + ieA,$2gP-T!^^^ + ]^e<^\^^-T'^^, (A12) 



From K^^,g — — F^^^^, and using eq. 



AlO 



we have 



2f 'V5 -I- A^g„.f"55 + A.f ^55 = 6$'g'^"T^,, (A13) 



13 



and hence, 



e* h^a 



^2f ^,„ + A^gA.f ^55 



■u5 



A r'' 

^i'^ -55 



(A14) 



These are 16 equations. Though the torsion in the Lh.s is antisymmetric, the r.h.s is a combination of symmetric and 



antisymmetric terms. Substituting for the 5D Levi-Civita connection from section (All, the above equation has the 
form, 



1 . cri/ ~ [C/o-Aiy — O^ActJ + [J(jA^ JiyAcrJ 

+ [e<p-^d5g,, + A^A^aggA. + A,A^95gA,] 



(A15) 



The first two terms in the r.h.s above are antisymmetric in a and v, while the third term is symmetric. The 
antisymmetry of torsion imphes that the symmetric terms in the r.h.s must be zero, and hence 



dsgai^ = 0. 



(A16) 



Consequently, 



rpO 



[d^A„ - d^A„] + [J^A^ - J^A^] . 



(A17) 



Thus the 20 equations of imposed condition (eq.fASl) have determined all 10 independent non- vanishing components 
of the torsion (equations All and A17|, and in addition imposed a constraint on the 10 independent components of 



the 4D metric gen, making them independent of x . 



The contorsion, connection and the Ricci tensor 



Some components of contorsion are directly prescribed by the imposed condition, namely K''-,- = K^, 



•5i 



pM 



i5- 



The remaining components of the contorsion can be calculated from the torsion components by using eq. [Tj Since 
T^j ■ are the only non vanishing components of torsion, it follows that 



K^.5 = K^5. 



• J5' 



K:i-, = -T:^ 



•i5 



i5 



g^^ (g,5T^5 + gssT^,,) 



K' 



5i 



K^5+T^5. 



K>:^,^-g^'{s.5T',x + Sx,T',.) 



K^,, = ^T:" 



:g'^(g.5T5,A + gA5T!^jv) 



(A18) 



Substituting for the nonvanishing torsion components given by equations All| and A17, along with the requirement 
that the 4D metric is independent of the fifth dimension (eq. A16|, the components of contorsion are found to be 






i5 









K. i5 + T. gj. 



K^,^ = e^^F'I^A,) + AaA, J'^e<&2 _ a^^^^^a^Aa^^), 



K^ 



-e$2A-F,(AA,) - (A,J")eA,AA$' + ;.(A,A" + e'^-^)d5ieA,Ax'S?^) - A,Ja + -F,a 



(A19) 



Now from cq. [5] we obtain all the connection coefficients. 



r. ^y — VpAy 



J /i-^Z^l 



■ 5^ 



C'sA/i + J^An, 



T^A T^A 

pM pM pM n 

^ -55 ~ ^ 1/5 ~ ^ -51^ ~ ^• 



(A20) 
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Taking i = ji and k — v va the Ricci tensor defined by eq. |15[ we have 



r 


- a.f-,^ 


+r:^.Ar".M 


_po- pA 


5m ■ 


- ^5? !^.^ 


+r^.Ar^5^ 


_f5 fA 
^•5A^i>M 






+r^.5r^.M 


"r^^gF. j,^ 






+f^5f^5M 


f5 f5 
~^-55^ -v^i 



(A21) 



Since f . is the same as V. when all the indices are four dimensional, the first line is clearly the 4D Ricci tensor. The 
terms in the subsequent lines can be re-expressed in terms of the 4D covariant derivative operator as follows 



a^fs 



n^ij — Hf^i, + V 1,1 . 5^ U^L . ^^ ± . ggX . j/^ T" J^ . 1/5^ . 5^ 



T^5 T^D 



r^ 



Substituting for the connection F. from eq. 



A21 



we find after some algebra this can be simplified as 
i?M^ = ^P^ + V^(55A^ + JgA^) - 55(V^A^ + J^A^) - J5(V^A^ + J^A^) + J^l^sA^ + JgA^) 

Similarly, the other components of the Ricci tensor are found to be 



-R^S = (^5^'^crti — 



i?55=0 



d^T^z = 



(A22) 



(A23) 

(A24) 
(A25) 
(A26) 



Note that neither the connection nor the Ricci tensor depend on the signature (e) of the fifth dimension. 
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Appendix B: Computing the modified Einstein tensor 

Here we provide sonic intermediate steps to go from equations |24| and |25| to eq. |26| and obtain a simplified expression 
for H^i,. With B^ = T^^^g, and H^,, defined by 



we show that 



J B. [g'^-^r:^^, - g'''''5r"„^] v^d'x = J R.jg'^-^ v^ 



H^, = V(^B,) - (V • B)g^, + J(^B,) - (J • B)g^,. 



d^x, 



Proof: 

Consider the first term in the integrand of the Lh.s of eq. IBI 



2g^"'^r!'^„ = g^^'-g-\d,5gxa + d^Sgx^ - dxSg^o.) + g''"id,gxc. + 5„gA^ - dxg^c.)Sg' 



u\ 



(Bl) 



(B2) 



(B3) 



The variations of the covariant metric in the above equation can be re-expressed in terms of the variations of the 
contravariant metric using the identity 5ga\ = — ga^gAi^'^g'^''- We note that 

g^<-g-^d^5g^^ = -g'^'d^gUh"'') - g^"5^g/3a(%'"') - ^^(^g''') 

g^"g^^5„JgA,, = g''"g'''a^<SgAa 

-g^^g^^SA^g^. = g'^^dxg.^Sg'^n + g''^dxg0.{5g^^) + g''^gp.dx{5g^P) 
Taken together, we obtain 

2g^"'5r''^„ = -2g''^9^gA.(<5g^^) - 2d^{5g^n + 2g'^^aAg^.(<5g'^^) + g''\p.dx{5g''^') - g^'^dxg^o.iSg"''). (B4) 



Next consider the second term in the integrand of the l.h.s of eq. Bl Since 2r"^ — g" S^gcA, 



25T\^ = (<5g"^)a^gaA +g"^9^('5goA) 
= -(<5g"^)a^g„A-gaA5^(^g"^) 
Using equations |B4| and |B5[ we find 

B. [g'-^SV^.^o. - g^''5T^o.^ ^/=g = B. [^d,{5gn + g'"^ga^5,,(<5g"^)] ./^ 

+ B^ 



(B5) 



-g^^a^gA.(5g'^'^) + Ig'^dxg.ASg'"') - \g^''dxg,a{5g^^) 



^ (B6) 



The first hne in the r.h.s of the above equation contains terms with the derivatives of the variation. We note that 
these terms are eventually going to be integrated. By integrating them by parts and ignoring the boundary terms, 
the above equation takes the form 



B. [g^'-SV^.o. - g'"''5r"aj y^ - [d,{B../^)5g^'' - 94g^''g„^B,y^)5g"^] 

+ B 



g^^5^gA.(<5g'^'^) + ^g'''9Agp.(^g"^) - ^g^'^5Ag^.(<5g^^^ 



To simplify the r.h.s of eq. |B7| it is useful to note the following identities 

^^(g'''') = g^'^g.A^MCg''^) = -g-''g^^a^(g.A). 



-g (B7) 



(B8) 



dug 



d,{V^) = ^x/^ = ^h^€'dM,) 



2g 



(B9) 



g'^a^(g,,) = g^^a^g^A) + g'^^a^(eA,AA$^) + 2g-^9^(eA,$^) + ga^5^(e$^) 
= g^^9^(g,A)+2J^ 



(BIO) 
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Using the above identities, the r.h.s of eq. |B7| becomes 

+ B.g„^g-''g^^9^(g.A)^g"^ - ^B,g'^-g„^g'^^9^(g^;,)5g"'^ 

- B,g''-9^(g,„)Jg"^ + iB,g--a,(g„^)Jg"'^ 
+ ^.i^Sg^" - B.J^g^'^g,;,^"^] y^. 
Rewriting the derivatives of the metric in terms of Levi-Civita connection, we find 

B. [g'^-'^r:^,, - g^'^^r-^^] V^ = %{^,)5g^^ ~ 9,(B,)g^^g„^<5g"^ 

- B,r!'^,jg'^^ + B.r^'^^g'^'^g^^.Sg"'^ 

+ B, J^Jg''- - B.J^g^^'g^^Jg"^] V^. 



(Bll) 



= [V^B, - (V • B)g^, + J,,B, - (J • B)g^,] ,5g'^''V=g (B12) 

Since the variation bg^^ is symmetric in the indices /i and v^ only the symmetric part of the r.h.s of the above 



equation will contribute to the equations of motion. Hence H^,^ will be given by eq. B2 

I 



